#This module contains the spatial, temporal, and soil hydraulic parameters
def spatialVariables_1D():
    totalDepth      = 0.32   # in meters ('m')
    axialNodes      = 17
    axialDistance   = totalDepth/(axialNodes-1)
    totalNodes      = axialNodes
    numOfActuators  = 1
    interPoints     = 500
    return totalDepth, axialNodes, axialDistance, totalNodes, numOfActuators, interPoints

def temporalVariable():
    samplingTime   = 1*6*4*60*60
    samplingTimeInternal = 1*6*4*60*60
    simulationTime = 1*6*4*60*60
    internalTimeSteps = int(samplingTime/samplingTimeInternal)
    totTimeSteps   = int(simulationTime/samplingTime)
    return samplingTime, samplingTimeInternal, internalTimeSteps, totTimeSteps


def Sandyclayloam():
    soilPars={}
    soilPars['thetaR'] = 0.10
    soilPars['thetaS'] = 0.39
    soilPars['alpha']  = 5.9
    soilPars['n']      = 1.48
    soilPars['m']      = 1-1/soilPars['n']
    soilPars['Ks']     = 0.000003639
    soilPars['neta']   = 0.5
    soilPars['Ss']     = 0.00001

    return soilPars